Universality of Velocity Gradients in Forced Burgers Turbulence 
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It is demonstrated that Burgers turbulence subject to large-scale white-noise-in-time random 
forcing has a universal power-law tail with exponent -7/2 in the probability density function of 
negative velocity gradients, as predicted by E, Khanin, Mazel and Sinai (1997, Phys. Rev. Lett. 78, 
1904). A particle and shock tracking numerical method gives about five decades of scaling. Using 
a Lagrangian approach, the -7/2 law is related to the shape of the unstable manifold associated to 
the global minimizer. 
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The universality of small-scale properties in fully de- 
veloped Navier-Stokes (NS) turbulence is frequently in- 
vestigated assuming that a steady state is maintained 
by a large-scale random force. For structure functions 
(moments of increments) universality with respect to the 
force is conjectured in the case of three-dimensional NS 
turbulence and proven for certain linear passive scalar 
models (see, e.g., Ref. [Q]). The universality of proba- 
bility density functions (p.d.f.) for velocity increments 
and gradients is a difficult question which, so far, has 
been mostly addressed within the framework of the pres- 
sureless model of Burgers turbulence, usually the one- 
dimensional Burgers equation 



dtu - 



(1) 



with white- noise-in-time forcing It is generally con- 
jectured that, when — > and the forcing is confined 
to large scales, the tail of the p.d.f. of velocity gradients 
^ at large negative values follows a universal power-law 
P{0 ^ The actual value of the exponent is how- 

ever a matter of controversy. Let us briefly recall some 
of the arguments found in the literature. 

A standard approach is based on studying the inviscid 
limit of the Fokker-Planck equation for the p.d.f. 

dtp - (fp) - + '^d^ m.m) P) - Bd^^p, (2) 

where the right-hand side expresses the diffusion of prob- 
ability due to the delta-correlation in time of the forcing. 
It was pointed out by Polyakov that the inviscid limit 
of (^) contains anomalies due to the singular behavior of 
the dissipative term vd(^ [{dxx£,\d P)- The value a = 3 is 
obtained if anomalies arc ignored y or if a piecewise lin- 
ear approximation is made for the solutions of the Burg- 
ers equation ||]. An operator product expansion (OPE) 
method borrowed from quantum field theory has been 
proposed for evaluating such anomalies and an argument 
presented in favor of a = 5/2 (actually, for velocity incre- 
ments and infinite systems) H. However, this expansion 
leads to a relation involving unknown coefficients which 



must be determined, e.g., from numerical simulations [3|, 
and restricts the possible values to 5/2 < a < 3 [7\. 
Anomalies cannot be understood without a complete de- 
scription of the singularities of the solutions, such as 
shocks, and of their statistical properties. For the case 
of a space-periodic system (as we shall assume), a cru- 
cial observation made in Ref. |^ is that large negative 
gradients stem mainly from preshocks, that is the cubic- 
root singularities in the velocity preceding the formation 
of shocks A simple argument was given in Ref. |^ 
for determining the fraction of space-time where the ve- 
locity gradient is less than some large negative value. 
This leads to a = 7/2 provided preshocks do not cluster. 
Determinations of the dissipative anomaly of (H) have 
been made by formal matched asymptotics [ p^ and by 
bounded variation calculus [0. With the assumption 
that shocks are born with vanishing amplitude from iso- 
lated preshocks, the value a — 7/2 was obtained 0, 0. 
Other attempts to derive a — 7/2 using also isolated 
preshocks have been made . Note that there are sim- 
pler instances, including time-periodic forcing |]l3j and 
decaying Burgers turbulence with smooth random initial 
conditions which fall in the universality class 

a = 7/2, as can be shown by systematic asymptotic ex- 
pansions using a Lagrangian approach. In the presence 
of forcing, the key issues which remained to be settled are 
the possible clustering of preshocks and, closely related 
to this, the possible birth of shocks with non-vanishing 
amplitude. The results presented hereafter almost com- 
pletely rule out such possibilities. 

Numerically solving the randomly forced Burgers equa- 
tion in the limit of vanishing viscosity in such a way as to 
obtain clean scaling for the p.d.f. of gradients represents 
a significant challenge. Broadly speaking, there are two 
classes of methods. On the one hand, methods involv- 
ing a small viscosity, either introduced explicitly (e.g. in 
a spectral calculation) or stemming from discretization 
(e.g. in a ffiiite difference calculation). Viscosity gives 
rise to a power-law range with exponent —1 at very large 
negative gradients [Q whose presence will make the in- 
viscid 1^1^" range appear shallower than it actually is, 
unless extremely high spatial resolution is used. On the 
other hand, there are methods which directly capture the 
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inviscid limit with the appropriate shock conditions such 
as the fast Legendre transform method of Ref. [|l5| (ex- 
tended to the forced case in Ref. This method is 
very well adapted to decaying Burgers turbulence with 
non-smooth Brownian-type initial data but, with 
spatially smooth forcing, it leads to delicate interpolation 
problems which have been overcome in the case of time- 
periodic forcing fl^ ; with white-noise-in-time forcing, it 
is difficult to prevent spurious accumulations of preshocks 
leading to a = 3. To avoid such pitfalls, we develop a La- 
grangian particle and shock tracking method [ p7| which is 
able to cleanly separate smooth parts of the solution and 
is particularly effective for identifying preshocks. The 
main idea of the method is to consider the evolution of a 
set of N massless point particles accelerated by a discrete- 
in-time approximation of the forcing with a uniform time 
step. When two of these particles intersect, they merge 
and create a new type of particle, a shock, characterized 
by its velocity (half sum of the right and left velocities of 
merging particles) and its amplitude. The particle-like 
shocks then evolve as ordinary particles, capture further 
intersecting particles and may merge with other shocks. 
In order not to run out of particles too quickly, the ini- 
tial small region where particles have the least chance of 
being subsequently captured is determined by localiza- 
tion of the global minimizer (see below). The calculation 
is then restarted from t = for the same realization of 
forcing but with a vastly increased number of particles 
in that region. This method gives complete control over 
shocks and preshocks 18 1 and allows an accurate deter- 
mination of the relevant statistical quantity while keeping 
a manageable number of degrees of freedom. 

Fig. |l| shows the p.d.f. of the velocity gradients in log- 
log coordinates at negative values for a Gaussian forc- 
ing restricted to the first three Fourier modes with equal 
variances such that the large-scale turnover time is or- 
der unity. Quantitative information about the value of 
the exponent is obtained by measuring the "local scaling 
exponent" , i.e. the logarithmic derivative of the p.d.f. 
calculated here using least-square fits on half-decades. 
It is seen that over about five decades, the local ex- 
ponent is within less than 1% of the value —7/2 pre- 
dicted by E et al. This value of the exponent was 
also obtained numerically (with a fewer particles) for 
other large-scale forcing instances with compactly sup- 
ported or exponentially decreasing spectra and also for 
non-Gaussian forcing (e.g. with Fourier amplitudes hav- 
ing a Bernoulli distribution or an uniform distribution in 
an interval). Evidence for non-clustering of preshocks 
is obtained by counting the average number of shock 
formations per unit time. For different types of large- 
scale forcing, we found that the typical mean number of 
preshocks per turnover time is comparable to the number 
of forcing Fourier modes significantly excited. For such 
forcings, the density of preshocks is found to vary by not 
more than 6% when the time step varies by two orders of 
magnitude around 6t = 10~^, which is hardly consistent 
with a power-law (and even a logarithmic) divergence as 



5t 0. Furthermore, we have checked that shocks are 
always born with vanishing amplitude (within numerical 
errors). 
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FIG. 1; P.d.f. of the velocity gradient at negative values in 
log-log coordinates obtained by averaging over 20 realizations 
and a time interval of 5 units of time (after relaxation of 
transients). The simulation involves up to A'^ = 10^ particles 
and the forcing is applied at discrete times separated by St = 
10~^. Upper inset: local scaling exponent. 

Turning now to theoretical results, let us briefly recall 
the construction of solutions developed by E et al. [^ , 
in terms of the dynamical system associated to the char- 
acteristics of (|^) in the inviscid limit . The force is as- 
sumed to derive from a Gaussian potential F{x,t), delta- 
correlated in time, periodic of period 1 and analytic in 
space. A statistically stationary regime is reached by 
taking the initial time at — oo. The central point of the 
construction is the following variational characterization 
of the solution at an arbitrary time (i = chosen for 
convenience) : 



d 

uix, 0) = — — min 
ox x(-) 



\x^{t)-F{X{t),t) 



dt 



, . (3) 

where the minimum is taken over all piecewise smooth 
(absolutely continuous) curves X{t) with t £ (— oo,0] 
such that X(0) = x. A curve minimizing the action 
m (|) is called a minimizer and should be understood 
as a fluid particle trajectory. It obviously has to satisfy 
for alH < the Euler-Lagrange equations: 



X{t) = C/(t), 
U{t) = f{X{t),t). 



(4) 
(5) 



Except for a flnite number of z-values, there exists a 
unique minimizer mM. The locations where there are 
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more than one minimizer correspond to shocks. The min- 
imizers converge exponentially fast backward in time to 
the trajectory of the unique fluid particle which is never 
absorbed by a shock. This trajectory is called the global 
minimizer because its action is minimal at any time; it 
corresponds to a hyperbolic trajectory of the dynamical 
system (^)-(|^). Associated to it, there are two curves 
in the phase-space {x,u): a stable (attracting) manifold 
r^*^' and an unstable (repulsive) manifold F^"). The min- 
imizers converge backward in time to the global mini- 
mizer and, thus, the graph of the solution is made of 
pieces of the unstable manifold with jumps at shocks. 
One of these shocks, called the main shock, is singled 
out. It is the unique shock which has always existed in 
the past (whereas generic shocks are born at some finite 
time t < 0); it may be shown that it corresponds to the 
position giving rise to the left-most and the right-most 
minimizers which approach the global one backward in 
time. The other shocks cut through the doublefold loops 
of the unstable manifold (see Fig. ||). We observe that 
their locations can be obtained by a Maxwell rule applied 
to those loops. Indeed, the difference of the two areas de- 
fined by cutting such a loop at some position x is equal 
to the difference of actions of the two minimizers defined 
by the upper and lower branches and, thus, vanishes at 
the shock location. 
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FIG. 2: Sketch of the unstable manifold F'"' in the (a;, u) 
plane at a time t with a preshock occurring. Shock locations 
are obtained by applying Maxwell rules to the loops. The 
velocity is shown as a bold line. 

We also observe that the structure just outhned has 
much in common with that appearing in the unforced 
Burgers equation. Indeed when f — 0, the solution to 
the Burgers equation can be constructed from the La- 
grangian manifold in the {x, u) plane, defined as the 
position and the velocity of fluid particles when ignor- 
ing shocks. This manifold is parameterized by the La- 
grangian coordinate a; denoting uq the initial velocity, 
we then have simply x = a + tuo{a) and u — uo{a). 
The actual solution with shocks is obtained by applying 
the standard Maxwell rule to the Lagrangian manifold. 
In the forced parameterization of the unstable 

manifold (e.g. by the arclength) is now the analog of the 
Lagrangian coordinate. But there are two important dif- 
ferences: first, in the unforced case, the time evolution 



of the Lagrangian manifold is explicit and linear while, 
in the presence of a force, the Euler-Lagrange equations 
(^)-(^) are not, in general, explicitly solvable and the un- 
stable manifold has a hyperbolic dynamic. Second, the 
smoothness of the Lagrangian manifold in the unforced 
case stems directly from the smoothness of the initial 
data, whereas in the forced case Pesin's theory must be 
used to show that when the force is indefinitely differen- 
tiable in space, so is the unstable manifold 

Using the smoothness of the unstable manifold, we 
now formally derive the —7/2 law, by an argument 
mostly borrowed from the unforced case 0. Let F(") = 
{{X{s), U{s))} with s real, be a parameterization of the 
unstable manifold at time t — 0. It is assumed for con- 
venience that s = corresponds to the global minimizer 
and that X'(0) > 0, where primes denote s-derivatives. 
The velocity is exactly obtained by eliminating from 
the unstable manifold the shaded areas determined by 
the Maxwell rules and the parts beyond the main shock 
(shown as dashed lines in Fig. ^). The surviving set of 
parameter values (excluding shocks) is denoted fi. Turn- 
ing to the statistical description, the p.d.f. of velocity 
gradients may be written 

p{O^{SidM^,0)-O)- (6) 

Because of homogeneity, we can integrate over the space 
period and then change from the x variable to the s vari- 
able to obtain 

P(0 = I' P{i)dx =l^jj{E_-^ x'ds'j . (7) 

Note that since a finite gradient is assumed, x cannot be 
at a shock position. Denoting by Sk the parameter values 
where the argument of the delta function vanishes, we 
obtain 

For very large negative values of ^, the s^'s must be near 
some s*j , corresponding to a local minimum of X'. Tay- 
lor expansions of X and U in the vicinities of the s*j's 
and the use of the Maxwell rule show that the s*j's are 
located in space-time near preshocks satisfying X' = 
and X" = with X'" > (see Fig. ||). Proceeding as in 
Ref. [0, we finally obtain, to leading order 

pio ^ e - (9) 

C = ¥ {ln.x"'>o \X'"\'/'\U'f/'6{X')S{X")dsyiQ) 

Hence, the constant involves the mean of 
at preshocks. Its evaluation requires the knowledge of 
the joint probability distribution of X', X", X'" and U' . 
From the Euler-Lagrange equations (^)-(^), we observe 
that a set of ordinary differential equations with non- 
linear stochastic forcing is easily obtained for X, U and 
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the aforementioned four variables. From these equations, 
using techniques similar to those developed in Ref. 
(where a subset of these stochastic equations is stud- 
ied), it should be possible, on the one hand, to make our 
derivation more rigorous (including for the non-clustering 
of preshocks) and, on the other hand, to obtain an up- 
per bound for the constant C in the —7/2 law. Note 
that the expression for C involves also an integral over 
the admissible set of parameters fl whose determination 
cannot in general be done by local analysis with ordinary 
differential equations. This is why only an upper bound 
is expected. 

As noted in Ref. the universality with respect to 
the forcing of the p.d.f. of large negative velocity gradi- 
ents may be extended to negative velocity increments, 
provided that they are not significantly influenced by 
shocks. Without understanding of all the mechanisms 
leading to small-amplitude shocks in the forced case, the 



issue of universality for the p.d.f. 's of velocity increments 
cannot be settled. A first step would be to determine 
numerically the distribution of shock amplitudes. Note 
that our technique may also be extended to the case of 
forcing at scales much smaller than the size of the sys- 
tem, a problem close to that considered by Polyakov [||, 
which is left for future work. 
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